Collective motions of a quantum gas confined in a harmonic trap 
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Single-component quantum gas confined in a harmonic potential, but otherwise isolated, is con- 
sidered. From the invariance of the system of the gas under a displacement-type transformation, it 
is shown that the center of mass oscillates along a classical trajectory of a harmonic oscillator. It is 
also shown that this harmonic motion of the center has, in fact, been implied by Kohn's theorem. 
If there is no interaction between the atoms of the gas, the system in a time-independent isotropic 
potential of frequency is invariant under a squeeze-type unitary transformation, which gives col- 
lective radial breathing motion of frequency 2vc to the gas. The amplitudes of the oscillating and 
breathing motions from the exact invariances could be arbitrarily large. For a Fermi system, ap- 
pearance of 2i'c mode of the large breathing motion indicates that there is no interaction between 
the atoms, except for a possible long-range interaction through the inverse-square-type potential. 

PACS numbers: 03.75.Ss, ll.30.Na, 03.75.Kk 



Experiments at ultracold temperature have stimulated 
theoretical efforts to explore the properties of quan- 
tum many-body systems in harmonic potentials. As 
the interaction between the atoms of the quantum gas 
is controllable 0, the density profile of A^-body non- 
interacting fermions has been investigated from the 
eigenstates of a harmonic oscillator For a one- 

dimensional harmonic oscillator, it is known that the 
time-dependent Schrodinger equation is invariant under 
the displacement- type transformation and squeeze- 
type transformation (up to a rescaling of time) . The 
invariance under the displacement-type transformation 
(DTT) gives wave functions whose probability distribu- 
tions move along classical solutions. The invariance un- 
der the squeeze-type transformation (STT) gives breath- 
ing motions to the distributions. For a time-independent 
potential of frequency t'c, it is numerically shown that 
the breathing motion has frequency 2i/c 

Harmonic motions of the centers of masses of the quan- 
tum gases have been clearly noticed in experiments and 
used to accurately measure the potentials (see, e.g., 0). 
Theoretically, dipole mode in the motion of the center 
has been found through various approximations, mainly 
based on the Gross-Pitaevskii (GP) mean-field formal- 
ism 1^ 0- In addition, for the time-independent po- 
tential, Kohn's theorem says the existence of exact ex- 
cited states with excitation energies 2'iTlhvc = 1, 2, • • • ) 
and further "harmonic-potential theorem" has been es- 
tablished 0. To my knowledge, however, the harmonic 
motion of the center has rather been understood from 
the insight that "collisions between atoms cannot alter 
center-of-mass momentum" |^ in the context of ultracold 
temperature physics. The implication of the breathing 
motion of the harmonic oscillator to the Bose-Einstein 
condensates has also long been studied in literature. In 
particular, it has been suggested that the GP equation of 
a time-dependent harmonic trap may be transformed to 
the equation of a time-independent potential through the 
unitary transformation of a harmonic oscillator and 
a scheme to understand the 2V(. mode is proposed [Tol |. 



Recently, a 2V(. mode has been found in an experiment 

El 

In this paper, single-component quantum gas confined 
in a time-dependent harmonic potential, but otherwise 
isolated, is considered, with an assumption that position- 
dependent part of the interaction between the atoms of 
the gas is written in terms of the differences of positions 
of two atoms. From the exact invariance of the gas sys- 
tem under a DTT, it is shown that the center of mass 
oscillates along a classical trajectory of a harmonic oscil- 
lator as it has been known for the time-independent po- 
tential through a different formalism |^ . It is also shown 
that this harmonic motion of the center has, in fact, been 
implied by Kohn's theorem. If the interaction between 
the atoms is absent, the system in a time- independent 
isotropic Vc potential is invariant under a STT up to 
a rescaling of time, which predicts the collective radial 
breathing motion of frequency 2vc of the gas. The am- 
plitudes of the oscillating and breathing modes from the 
exact invariances could be arbitrarily large, as they are 
determined purely by the classical solutions. In this re- 
spect, the exact modes may be different from other modes 
found in the linear approximations 0- For a Fermi sys- 
tem, appearance of 2i'c mode of the large breathing mo- 
tion signals that the gas is in the region of no interaction 
between the atoms, except for a possible long-range in- 
teraction through the inverse-square-type potential. The 
collective harmonic motions of the fermions will show 
that the complete set of the numbered coherent states 
[T^ 0, 01 of a harmonic oscillator can be used in stack- 
ing the fermions, with equal validity, as the eigenstates 
are used. For simplicity, the nonlinear Schrodinger (gen- 
eralized GP) equation will be considered first, and the 
formalism will be extended to the many-body systems. 

A D-dimensional nonlinear Schrodinger equation 
(NLSE) with a time-dependent harmonic potential is 
given by 



0(t,u;(t))vI/(r,t)+.g|vl/(r,t)|2«vI/(,-^,t) ^0, (1) 
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where V — {-^r—, ir—-, • • • ir—)-, fn is the mass of an atom, 
and (7,n are constants. When fi satisfies the Hill's equa- 
tion, the classical equation of motion of a harmonic os- 
cillator. 
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with the overdot denoting differentiation with respect 
to t, the displacement-type unitary operator for the D- 
dimensional harmonic oscillator can be given as 
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where 5i is defined through the relation 

Making use of the operator, one can find that 
^(r,t) = Uf{r,t)'^{r,t), 



(4) 
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^{r-f,t), (6) 



with / = (/i, /2, • ■ ■ I /d)- If the unitary operator Uf{r, t) 
is applied on (1), from the relation 

Uf{r,t)0{t,w{t))u]{r,t) = 0[t,w{t)) (7) 

and the fact that [//(r, t) [|*(f, i)|2*(r, i)] = 

one can find that the transformed 
equation is the same one with the replacement of 'I'(r, t) 
by ^(r, <). If ^'(r, t) is a solution of (1), vE'(r, i) is also a 
solution of the same equation, and the number density 
|\E'(r, i)P of the transformed solution has harmonic 
motions with respect to the original one, without 
changing the shape. If Wi{t) is a constant Wd and 
the transformation is applied on a stationary solution 
of (1), then the number density of the transformed 
solution oscillates sinusoidally along the j-th direction. 
For D — 1, the invariance has been known through a 
different formalism 

The invariance under the DTT can be easily ex- 
tended for a A^-body system, if the interaction 
F(r(l),r(2), • • • ,r{N)) between the atoms satisfies 



y(f(l) + c, f(2) + c, • • • , f{N) + c) 
= V{f{l),r{2),--- ,f{N)), 



(8) 



where f(j) denotes the position of the j-th atom and c is 
a constant vector. By defining Un as 
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one may find the relation 

UNON{t,w{t))Ul = ON{t,w{t)), 



(10) 



where 



ON{t,w{t)) EE -ih— + HN 
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For a given time, the operations of Un and Uf{'r,t) 
amount to moving the space coordinate along the vector 
/ to find new solutions, up to the multiplicative phase 
factors. The probability distribution of the transformed 
wave function of the A^-body system will thus have the 
harmonic motion with respect to the distribution of the 
original wave function, as the number density of the 
transformed solution of the NLSE does. (10) is valid for 
both of the Fermi-Dirac and Bose-Einstein statistics, and 
is an exact explanation for the longstanding observation 
that the center of mass of the single-component gas con- 
fined in a harmonic potential moves harmonically, while 
the harmonic motion of the center has been known for 
the time- independent potential Q. 

For the case of constant frequency Wi{t) = Wd (i = 
1,2,--- ,-D), one may find that the harmonic motions 
of the centers of masses have, in fact, been implied by 
Kohn's theorem For this, by letting 
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with real constants zoi, (pi, one can find the relations 
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e ^ exp 



= e**' exp[z,(t)al(j)-<W«.(j)], 

where 

z,it) = zo^e-'^^'"'+'^'\ 



(13) 
(14) 



d 



mwci dxi{j) 




and z*{t) denotes complex conjugate oi Zi{t). In (15,16), 
(pi is a real constant coming from that 5i is defined up 
to a constant, and from now on we will set (pi — Q {i — 
1, 2, • • • , £)). One can thus find that J7jv is written as 



Un = \{e^v(z,{t)A\ - z*{t)A,^ (17) 
i=i 

D 

= ne-^^o./2exp(z,(t)4)exp(-z*(i)A0, (18) 
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where 



N 



A, = ^ai(j)• 



(19) 



Indeed, in Ref. it has been proven that 

If |G > is the many-body ground state of with energy 
eigenvalue Eg, a wave function of the system is given as 



D 

^ e 

i=l 



e~TiEGtJJ^-Nzl/2 



\G > 



(20) 



while, as suggested by Kohn's theorem ^^1^ jG > is 
an eigenstate of the Hn with the eigenvalue Eq + Ihwd 

From now on, isotropic potentials will only be consid- 
ered, so that Wi{t) — w{t) for all i. For the STT, the 
unitary operator 



n 



mwcrfit) 
X exp 



exp 



i mr]{t) 
2h r]{t) 



1 ( mwc7f{t) ^ 



(21) 



is introduced, with positive constants il and Wc {— 2t:i>c)- 
For the time being, it will be assumed that ri(t) is an 
arbitrary smooth positive function of t. If Us is applied 
on x(r, T(t)), the transformed function will be 



*(r,i) = Usxir,rit)) 

n 



XX( 



exp 



mf][t) 



2h r]{t) 



n f 



mwc ri{t) 



.r{t)). 



(22) 



A NLSE with a time-independent harmonic potential is 
given as 



OMt),w,)xir,rit)) + g\x{f,Tm^"xir,rit)) = 
where 



OciT{t),Wc) 

From the fact 



-ih 



dr 2m 



(23) 
(24) 



Us{\x{r,r)\^''xir,r)) 

nD 

n 



mWcTf' [t) 



|*(f,i)|^"vl/(f,t), (25) 



one may find that the NLSE of (23) can be transformed, 
through the STT, to that of (1) only when 



(26) 



(27) 



dt mWcTj'^ 

Further, if ri{t) satisfies 

mi] + mw (t)r] — 0, 

mrj"^ 

one can find the relation 

UsOAT{t),w,)Ut - (^^) Oit,wit)) (28) 

which is well-known for the one-dimensional harmonic 
oscillators 0, ^3- When r satisfies (26), by applying 
the operator Us on (23), one can find that the NLSE is 
transformed as 

Us [OMt),w,)x{r,T{t))+g\x{r,T{m'-x{r,r{t))] 
mWcTj^it) 



1-^ 
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mWcr]-^ (t) 



|*(f,i)|2"*(f,<) 



0. (29) 



Through the STT, (23) is thus transformed to (1) when 

nD ^ 2. (30) 

As far as a NLSE is concerned, this condition is identical 
to that given in Ref. found in a different context 

using a trial function. For n = 1,2, the squeeze-type 
relation (29) has been suggested in Ref. 0. 

In the harmonic oscillator, it is known that r]{t) and a 
constant fl are written as (see, for example, Refs. 0,0|) 

r/(t) ^ ^/u^{t) +v^{t), n^m[v{t)u{t)-u{t)v{t)], (31) 

with two linearly independent real solutions u{t) and v{t) 
of 



Xcl + w'^{t)Xcl = 0. 



(32) 



From now on, we only consider r(i) and r]{t) satisfying 
(26,27). When (30) is satisfied, the unitary relations can 
be used to find a solution ^'(r, i) of (1) from a known 
solution xii^jTit)) of (23), as 

^{r, t) = Uf{r, t)Usx{r, T(i)). (33) 

If e^*^'^''*'xo(^) is a solution of (23) with a constant /i, a 
solution of (1) is given as 



*(f,i) 



mwc'rf'{t) 

\i=i 



u{t) — iv{t) 



-{5i + mfiXi) 



X exp 



imfj{t) 
2hri{t) 



(r- I f 



Xo( 



» r-f 
mWc rj{t) 



). (34) 
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For the case of / = 0, one finds that 

/ V D 



mwc 'ri{t) 



lxo( 



mWc Tj{t) 



)?■ (35) 



For a given time with f — 0, (35) shows that the number 
density of ^'(r, i) is given from that of e^'^'^(*)xo('') by 
globally rescaling the space coordinate along the radial 
direction, with a multiplication factor needed to keep the 
total number of particles. 

The squeeze-type unitary relation can be extended for 
the interacting A^-body system, if interaction between the 
atoms is of homogeneous degree —2, so that 



V{ar{l), af(2), • • • , af{N)) = a-^V{f{l)f{2), 
with a constant a. By defining 



,r{N)) 
(36) 
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2m 
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X exp 
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II exp 



i mf]{t) _,2 



2h ■q{t) 
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(38) 



one may find the relation 

UN,sON,c{r{t),Wc)U],^^ = ON{t, w{t)). (39) 

If XN{T{t)) [= e-^^-(*)/^Xiv,o(r(l), r(2), • • • , r{N))] sat- 
isfies the Schrodinger equation 



OM,rMt),Wc)XN{T{t))=0, 

a solution of the Schrodinger equation 

OM{t,w{t))-9N{t)^0 

is thus given as 
*w(t) 

u{t) — iv(t) 

~~W) 



(40) 



(41) 




E/h 



imr]{t) 



2hri{t) 

xxw,o(r(l),f(2),--- ,f(7V)), 



where f(j) 



n r(j) 



As in the number density of the 



NLSE, for a given time, the probability distribution of 
\I'jv(i) is found from that oixN{T{t)) by globally rescaling 
the space coordinate along the radial direction with a 
multiplication factor. 

When w{t) = Wc (and (30) is satisfied for the NLSE), 
the squeeze-type operators in (21,38) transform (23,40) 
into themselves with a replacement of r by t (up to a 
rescaling of the time) , respectively, and thus the unitary 
relations give invariance. If ri{t) is denoted as ric{t) in 
this case, r]c{t) is written without losing generality as 



Vc{t) = 



mWc 



(42) 



X ^Ja^ cos^ wc{t - to) + A-^sin^ Wc{t - to), (43) 

with a constant to and a non-zero constant A. ric{t) is 
a periodic function of time with frequency 2i/c, and thus 
the probability distribution of the wave function obtained 
from a stationary state of the A^-body system through the 
STT has the breathing motion of frequency 2i'c- 

The invariance under the STT is possible only when 
the atoms interact each other through the potential sat- 
isfyin g ( 27) . Such potentials are the inverse-square- 
type [13 and the two-dimensional Dirac delta poten- 
tials, including other variants from the one-dimensional 
Dirac delta potential. Since Dirac delta potential can- 
not be used for the interaction of fermions due to 
the exclusion principle, the Fermi system with the col- 
lective 2i'c breathing mode should be interpreted as 
that of no interaction between the atoms (except for 
a possible long-range interaction through the inverse- 
square- type potential). In the Tonks-Girardeau limit 
of the Tonks gas, interacting bosons behaves like non- 
interacting quasifermions, and the mean-field description 
has been given by Kolomeisky et al. [T^ . Since non- 
interacting particle systems in one dimension have the 
invariance, the mean-field equation should also be invari- 
ant under the STT, and the equation of Kolomeisky et 
al. satisfies the requirement of (30). 

In order to realize the large breathing mode from a 
stationary state in a time-independent trap of frequency 
Vc, one may modulate w{t) until r]{t) determined by (27) 
has a large oscillating behavior. In this case a convenient 
choice for the arbitrary constant is J7 = mWc, and the 
limiting case of v4 — > or A ^ 00 of (43) shows the exis- 
tences of such modulations. The modulation of frequency 
has already been widely used in experiments (see, e.g., 
[lll|). If the amplitude of the center-of-mass oscillation 
decays in an experiment of a single-component gas using 
a time-independent harmonic potential, it indicates that 
the gas is not completely isolated. While a 2i'c breathing 
mode is found and a sinusoidal curve is used to fit the 
size of the breathing gas in Ref. [H, (35,43) imply that 
the difference between the radius of the breathing motion 
from the invariance and the sinusoidal curve will be clear 
for j4 >> 1 (or A << 1) aroimd the time the gases are 
most compressed. 
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In summary, isolated single-component quantum gas in 
a harmonic potential is considered. The invariance under 
the DTT shows that the center of mass of the gas moves 
along a classical trajectory of a harmonic oscillator. For 
a time-independent potential, it is shown that the har- 
monic motion of the center has, in fact, been implied by 
Kohn's theorem |^ 0|. For a Fermi system, through 
the invariance under the STT, it is shown that appear- 
ance of the large radial 2vc breathing mode indicates that 
the atoms do not interact with each other (except for 
a possible long-range interaction). While some NLSEs 
and noninteracting Fermi gas in any dimension are in- 
variant under the STT, the collective breathing motion 
of the one-dimensional gas can be understood from the 



evolution of an eigenstate of the harmonic oscillator into 
a coherent state when modulation is applied to the fre- 
quency. Alternatively, in one dimension, the collective 
breathing and oscillating motions in a time-independent 
potential can be interpreted as the result of using gener- 
alized coherent states [ij, Hj, 13 of a harmonic oscillator 
in stacking fermions. 
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